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Starting from a microscopic model, we investigate the optical spectra of molecules in strongly- 
coupled organic microcavities examining how they might self-consistently adapt their coupling to 
light. We consider both rotational and vibrational degrees of freedom, focusing on features which 
can be seen in the peak in the center of the spectrum at the bare excitonic frequency. In both 
cases we find that the matter-light coupling can lead to a self-consistent change of the molecular 
states, with consequent temperature-dependent signatures in the absorption spectrum. However, 
for typical parameters, these effects are much too weak to explain recent measurements. We show 
that another mechanism which naturally arises from our model of vibrationally dressed polaritons 
has the right magnitude and temperature dependence to be at the origin of the observed data. 


I. INTRODUCTION 

When matter is strongly coupled to light, the interac¬ 
tion cannot simply be thought of in terms of absorption 
and emission processes. Instead we must consider the 
eigenstates of the fully coupled matter-light system. The 
paradigmatic example of this is the existence of exciton 
polaritons, hybrid matter-light particles formed by the 
strong interaction between excitons and photons P, [2| . 
Matter-light coupling can be engineered by confining 
light in optical cavities, so as to modify the density of 
states and the coupling to matter. For weak coupling, or 
a bad cavity, cavity losses are fast so one can eliminate 
virtual processes where photons are in the cavity. This 
gives Fermi’s golden rule, but with the cavity density of 
states modifying the emission rate, as first discussed by 
Purcell Q. When coupling is strong, first order pertur¬ 
bation theory (i.e. Fermi’s golden rule) fails, as there can 
instead be coherent emission and reabsorption of photons 
before light leaks out of the cavity ii- 

A natural context in which strong matter-light cou¬ 
pling arises is between organic molecules and light in 
semiconductor inicrocavities. Because of the existence of 
conjugated tt bonds in organic molecules, electronic tran¬ 
sitions can acquire large dipole moments , leading to 
very strong coupling to light. When such molecules are 
placed in optical microcavities this leads to huge polari- 
ton splittings [Ml. These scales allow such experiments 
to be performed at room temperature, whereas for many 
inorganic materials, cryogenic temperatures are required. 
The polariton splitting is due to a collective phenomenon: 
the electronic transitions of many molecules couple to 
radiation, and as such the polariton splitting grows as 
the square root of the molecule density. In contrast, in 
weak coupling, the rate at which one molecule emits is 
independent of whether or not any other molecules are 
present (T^ . 

Much of the recent work on organic microcavity po¬ 
laritons (see, e.g. Ref. [l^ for a recent review) has been 
focused on condensation and lasing (T5l - [l8j . involving a 
strongly pumped system, and the appearance of macro¬ 


scopic quantum coherence. There has however also been 
significant recent work on the effects of matter-light cou¬ 
pling in the vacuum state, i.e. without strong pump¬ 
ing. Such work aims to understand how the physical 
and chemical properties of organic molecules are affected 
by strong coupling to electromagnetic modes. Examples 
of this include modifying the rates of photochemical re¬ 
actions [l^ , or modify ing the transport properties of or¬ 
ganic semiconductors More recently, there has 

also been experimental and theoretical [241 - 1^ work 
on coupling the vibrational state of organic molecules 
to infra-red radiation, leading to molecular optomechan- 
ics. Theoretical work has also studied how strong 
matter-light coupling to electronic states can suppress 
the effects of disorder and vibronic features in the po¬ 
lariton spectrum. Of particular interest for the pre sent 
paper is a recent work from the Ebbesen group (28[, in 
which the optical spectra of strongly-coupled organic mi¬ 
crocavities were studied by varying molecular concentra¬ 
tion and temperature, and paying particular attention to 
the relative weights of the resonant features in the ab¬ 
sorption spectra: the two polarit on p eaks, and a third 
peak at the bare excitonic energy |^ . 

Our aim in this manuscript is to examine the behavior 
of such strongly-coupled organic microcavities starting 
from various microscopic models, allowing quantitative 
predictions of the extent to which a self-consistent adap¬ 
tation of the molecular state, driven by coupling with 
light, may occur. 

To understand the variation of the optical spectra with 
both concentration and temperature the models which 
we consider all contain a variable degree of coupling to 
light. This is because a molecule that has strictly zero 
coupling to light is not visible in the absorption spectrum, 
while molecules with a small but non-zero coupling will 
lead to absorption at the bare molecular energy. In or¬ 
der that the coupling to light can vary self consistently 
(in response to the Rabi splitting), it must depend on 
some adaptable feature of the state or environment of 
the molecule, i.e. there must be some physical property 
that can vary, which determines the strength of matter- 
light coupling. We refer to this concept hereafter as 
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“self-consistent molecular adaptation”. We refer to this 
process as “self-consistent” because the effective matter- 
light coupling depends on (some aspect of) the molecular 
state, and the molecular state is modified because of how 
its energy depends on the matter-light coupling. In the 
first part of this manuscript we investigate in detail two 
candidates that could lead to self-consistent adaptation: 
rotational and vibrational degrees of freedom, we also 
consider an extension of these models to generic (classi¬ 
cal) aspects of the molecules physical or chemical state. 
While we do find a temperature dependence of the op¬ 
tical spectra, the involved energy scales turn out to be 
incompatible with the observation of Ref. [2^ . In the fi¬ 
nal part of the present work, we examine how our model 
of vibrationally dressed polaritons naturally predicts an 
effect whose energy scale is of the right magnitude to 
explain the data. This effect does not involve the renor¬ 
malization of the coupling strength, but instead involves 
the effect of vibrational replicas and their coupling to 
the excitonic transition on the optical spectra. Our re¬ 
sults thus show a first example of the rich, and presently 
poorly understood behavior that can stem from the inter¬ 
play of strong matter-light coupling with strong coupling 
to vibrational/conformational modes of the molecules. 

We start by noting that the existence of a peak at 
the bare energy of the exciton, brought forward as ev¬ 
idence of novel physics in Ref. [l^, is not unexpected; 
such a “residual excitonic peak” has been seen in many 
cases, for example Ref. discussed theoretically, and 
demonstrated experimentally, the appearance of such a 
feature in a GaAs/AlGaAs heterostructure, containing 
quantum wells inside a DBR microcavity. While such a 
peak comes from the spectral weight of the exciton line, 
it is important to note that this peak cannot be viewed 
simply as excitons which do not couple to light: if they 
did not couple, they would not be visible in the absorp¬ 
tion or transmission spectrum. The origin of the peak 
can be understood physically as coming from the subra¬ 
diant excitonic states due to inhomogeneous broadening. 
In a disordered system, the coupling to the photon mode 
picks out a specific superradiant state, which forms the 
polaritons, while the other states — orthogonal to this 
superradiant state — remain at the bare exciton energies. 
However, because of energetic disorder, the superradiant 
state is not an energy eigenstate, and a residual coupling 
between the superradiant and subradiant states exists, so 
that the spectral weight of the subradiant states is visible 
in the optical spectrum 

A simple analytical treatment shows that the weight of 
this residual excitonic peak does decrease as the matter- 
light coupling increases. The existence of this peak is 
thus consistent with the behavior seen in Canaguier- 
Durand et al. [2^. However, on its own this explana¬ 
tion cannot account for the temperature dependence ob¬ 
served, as the residual excitonic peak should be unaf¬ 
fected by temperature, unless ksT approaches optical 
energies, of the order of leV (for comparison SOOKr; 
25meV). One of the main goals of this paper is to ad¬ 


dress how this temperature dependence may occur. 

As will become clear in the following, to discuss the 
microscopic theory of such effects, it will be crucial 
to consider the physics of ultrastrong matter-light cou¬ 
pling [H, 113 , and the breakdown of the rotating wave 
approximation (RWA). This requires retaining “counter¬ 
rotating” terms in the matter-light coupling Hamilto¬ 
nian. These terms, which involve simultaneous creation 
of pairs of excitations, are typically considered to be non- 
resonant and so are often neglected. However, if the 
matter-light coupling is a signihcant fraction of the bare 
exciton and photon energies, then these terms have a 
non-negligible impact. Such behavior has been seen in 
both inorganic [13 and organic [Ill[ll,ll3| systems, with 
a current record of a coupling strength 87% of the bare os¬ 
cillator frequency [38l |- Our focus in this paper is on the 
more typical regime where such counter-rotating terms 
cannot be neglected, but remain sufficiently small to be 
treated perturbatively. 

The rest of the paper is structured in two main sec¬ 
tions. In section El we consider how temperature de¬ 
pendence can arise due to self-consistent adaptation of 
the rotational and vibrational degrees of freedom of the 
molecules, via a mechanism very similar to that pro¬ 
posed in Ref. [13 . For the orientational degree of free¬ 
dom, we consider both free molecules, and molecules with 
randomly pinned orientations as appropriate in a poly¬ 
mer matrix. We will see that in such systems we do 
predict a temperature dependence of the residual exci¬ 
tonic peak. However, while this effect could potentially 
be observed in other experimental realizations, the en¬ 
ergy scales (temperatures) required and the scaling with 
molecular concentration are not compatible with the ex¬ 
perimental observations reported in Ref. [13 • 

In section El we instead consider a different effect, 
arising from the interplay of vibrational modes with the 
matter-light coupling, which is able to reproduce similar 
behavior to that observed in experiments. Specifically we 
find that vibrational excitations dress the residual exci¬ 
tonic peak in a strongly temperature dependent manner. 
Moreover, the form of the vibrational dressed spectrum 
shows that the spectral feature at the exciton energy can 
have a more complex interpretation than that previously 
considered [13 ■ 

A brief but self-contained account of the main theoret¬ 
ical methods used throughout this paper is given in the 
appendices. 


II. SELF-CONSISTENT MOLECULAR 
ADAPTATION DUE TO STRONG COUPLING 

In this section we consider whether self-consistent 
molecular adaptation can enhance matter-light coupling 
by renormalizing the bare matter-light coupling strength. 
We consider models in which the effective matter-light 
coupling strength of a given molecule depends on the 
configuration of that molecule, such as its orientation, or 
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its vibrational state. We then ask how this same matter- 
light coupling modifies the energy landscape for the aux¬ 
iliary parameters describing the configuration. This leads 
to the idea of self consistency — if strong coupling leads 
to a reduction of the ground state energy, the energy 
landscape is deformed so as to favor auxiliary param¬ 
eters for which the effective matter-light coupling is as 
large as possible. Our aim is to derive this from a micro¬ 
scopic model, and so quantify this effect. In the following 
we consider two potential scenarios involving adaptation 
of either orientational or vibrational degrees of freedom. 

If such a self-consistent enhancement of matter-light 
coupling occurs, then this can lead to a temperature 
dependent effective coupling, and thus to a tempera¬ 
ture dependence of the residual excitonic peak. We 
show that such an effect exists, but that its strength 
is relatively weak, and that the relevant energy scale 
shows no collective enhancement, i.e. the presence of Nm 
molecules does not lead to a Nm enhancement of this en¬ 
ergy scale, because it must compete with the extensive 
entropy gain from orientational or vibrational disorder. 
As such, while increasing the molecular concentration 
will increase the polariton splitting, it has little effect on 
the self-consistent orientation. Changing the bare oscil¬ 
lator strength of the molecules does however affect both 
the polariton splitting and the self-consistent molecular 
adaptation energy scale. 

Before introducing any auxiliary variables, the basic 
Hamiltonian which we consider is an extended Dicke 
model, including diamagnetic terms: 


Cn 2 : 
2 


k n 

+ Xl5k,n<Pk(rnX + ( \/^‘/5k(r„) 


( 1 ) 


where the field :/Jk(r) = -I- is written 

in terms of the bosonic creation and annihilation oper¬ 
ators describing photon modes labeled by their 

in-plane momentum k, and energy Wk- The Pauli matri¬ 
ces describe the electronic state of the molecules. 

For completeness, we therefore also included the diamag¬ 
netic terms, arising from the term of the minimal 
coupling Hamiltonian [d^. 

Furthermore, in order to consider varying the cavity 
mode volume while respecting the Thomas-Reiche-Kuhn 
sum rule we use t/k.n = y/^nD^fn where 0 < /„ < 
1 is the oscillator strength of the given molecule. The 
coefficient Uk depends on the electric field strength of a 
single photon, and the properties of the effective charges 
that respond to the field. 

If we assume an uniform distribution of molecules, mo¬ 
mentum conservation can be used to write the diamag¬ 
netic term in Eq. o as NmJ2k 

This term can then be removed by a Bogoliubov transfor¬ 
mation —>■ cosh^ktf’k+sinhdk'^ki yielding the effective 


Hamiltonian: 


= ^WkV'k^k+ (2) 

k n 

where the on-site Hamiltonian is given by 

hn = y < + 9k,n<p{rn)(J^, (3) 

k 

and the renormalized parameters Wk,3k are: 


^k — \/a^k(a^k T dAf^T^ZIk); 



(4) 


with Nm the number of molecules in the mode volume. 

In the dipole approximation, for a molecule with a sin¬ 
gle mobile electron Uk = C/^kV with V the mode vol¬ 
ume and C = Tt?/{A nrireo) quantifying the electronic 
response of a single electron in terms of the vacuum per¬ 
mittivity So and its reduced mass m^. For molecules 
involving many conjugated bonds, the coefficient C is re¬ 
placed by a sum over all mobile charges. It is important 
to note that changing the cavity length changes both the 
mode volume V (and hence both T>k and 5k,n) and the 
spectrum of photon modes Wk- In the following numeri¬ 
cal results we will fix the values of the polariton splitting 
gk,n^/Nm, exciton energy e„, and coupling strength /, 
and use these to determine D\^Nm- 

Before considering the role of vibrational and orien¬ 
tational degrees of freedom, we consider how the exis¬ 
tence of a temperature dependent matter-light coupling 
strength, 5eff,k,n would be seen in the absorption spec¬ 
trum. For this purpose, it is sufficient to consider the 
absorption spectrum of Eq. and its dependence on 
5k,n. Appendix lAl summarizes how the absorption spec¬ 
trum can be calculated, including the counter-rotating 
terms in the matter-light coupling (^ . In performing 
these calculations, as noted above, it is necessary to in¬ 
clude disorder in order to see the residual excitonic fea¬ 
ture. We will consider disorder in the exciton energies, 
denoted by the energy distribution h(e) Eor simplicity we 
consider only disorder in the energies and we ignore the 
subleading effect of the exciton energy distribution on 
the coupling strength gk,m by using the averaged value 
5 ^ = J2n9knl^m- As discussed in the appendix, we 
consider the quantity ak(j^) = in terms 

of the retarded Green’s which is proportional to 

the absorption spectrum for a good cavity. The Green’s 
function has the form: 

Gkxxi^) = ^ ! \ - ■- -/ X ; (5) 

V — -\- 2u}k^k,xx{l^) + 

where the excitonic self energy for Eq. ©, can be written 
as Ek,a:x = -I- with: 


^RWA 
^k, + - 


(^) 


-9lNr. 


/ oo 

deh{e 

-OO 


tanh(/3e/2) 
V + 17 + — e 


( 6 ) 
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Here /3 = {kBT)~^ is the inverse temperature and 7 is the 
homogeneous linewidth of the excitons. In the following 
we will present results both for 7 = 0'*' and for small but 
non-zero 7 as indicated in the figure captions. 

The spectrum is shown in Fig.[Tl focusing on the resid¬ 
ual excitonic peak, to show its dependence on the polari- 
ton splitting. As noted above, such a feature has been ob¬ 
served and commented on several times before, e.g. [ 13 - 
[ 3 ^. When the polariton splitting g\/Nm is increased, 
the splitting between the lower and upper polaritons in¬ 
creases, and the exciton spectral weight of the feature at 
the exciton energy decreases. The asymmetry between 
the shift of the lower and upper polaritons arises due to 
the diamagnetic term Dk renormalizing the photon en¬ 
ergy. If the coupling to light is weak, so that the RWA 
is valid, this asymmetry vanishes. In Fig. [T] (b) we show 
how the spectrum is modified by including the effects 
of cavity losses and non-radiative excitonic decay as de¬ 
scribed in Appendix The main effect these processes 
have is to broaden and thus reduce the height of the po¬ 
lariton peaks, there is also additional broadening to the 
central peak. 

Because the only temperature dependence of Eq. (EH) 
is via the combination tanh(/Ie/2), the spectrum is tem¬ 
perature independent while ksT <C e (with e of the order 
of the exciton bare energy). Therefore, as noted in the 
introduction, the experimentally observed temperature 
dependence cannot occur from this mechanism alone, un¬ 
less the effective value of gk,n is made temperature de¬ 
pendent via its dependence on configuration. We now go 
on to consider if this effect is plausible. 

To consider molecular adaptation, the Hamiltonian in 
Eq. (H will be modified to include either orientational or 
vibrational degrees of freedom. These are illustrated in 
Fig. [21 We next introduce these modifications and then 
discuss how they may be treated. 

The Hamiltonian including an orientational degree of 
freedom takes the same form as Eq. (H with the on-site 
part 

hn — ^niPn) T ^ ^ ffk.n COs(^yt)(^k ; (^) 

k 

where the angle 9n parametrizes the orientation of the 
dipole moment of the nth molecule with respect to the 
polarization of the cavity electric field, thus reducing 
the oscillator strength. The term A„(0„) represents the 
bare dependence of the Hamiltonian on orientation. This 
term allows one to model pinning of the orientation 0„. 
Eor simplicity we consider only a classical orientational 
degree of freedom; the corresponding quantum theory 
would require us to also include a rotational kinetic en¬ 
ergy term, and diagonalize the resulting Hamiltonian. 
The effective matter-light coupling strength, which de¬ 
pends on the distribution of angles On adopted by the 
molecule, can be written as „ .eff = -glnii^^OsHOn))). 
where the double angle brackets represent both an en¬ 
semble and thermal average. 



V [eV] 

FIG. 1. (color online) Evolution of k = 0 absorption spec¬ 
trum, a{u) vs polariton splitting gk=oy/Nm for the system 
with no auxiliary degrees of freedom. Plotted for ksT — 
0.025eV (i.e. T = 290K), / = 0.05,aJk=o = 2.1eV, and a 
truncated Gaussian distribution h{e) oc 0(e)e^'“'^°^ with 
£0 = 2.0eV, a = O.OleV. Panel (a) shows the results with a 
perfect cavity while (b) includes the effects of cavity losses at 
rate k = 0.075eV and excitonic non-radiative decay at rate 
7 = 10~^eV also the width of the energy distribution is larger, 
a = 0.025eV as discussed in Appendix m 

To consider vibrational degrees of freedom, we again 
start from the transformed Hamiltonian, Eq. (ED, and 
now consider the following modification: 

hn = y < + ^ 5k.nV5k(r„X-b 

k 

+ • ( 8 ) 

m ^ ' 

Here bn^m,bh m describe the mth harmonic vibrational 
mode of molecule n, with the mode having frequency 
flm and its coupling to the electronic state being 
parametrized by the Huang-Rhys parameter y/Sm- In 
this case, defining the effective oscillator strength is more 
involved: the effective oscillator strength depends on the 
matrix element describing the overlap between the vibra¬ 
tional states in the ground and excited state manifold. 
We return to this point in later sections. 

Eor both the orientational and vibrational degrees of 
freedom, our aim is to find how the matter-light cou¬ 
pling is self-consistently modified by these auxiliary de¬ 
grees of freedom: i.e. how the presence of matter-light 
coupling modifies the distribution of orientational or vi¬ 
brational states, and how this in turn affects the effective 
matter-light coupling strength. We consider a case with¬ 
out any strong pumping, and with a temperature such 
that ksT <C w, e, which is typically satisfied even at room 
temperature for organic polaritons. As such, the origin 
of the “self-consistent” dependence of configuration on 
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FIG. 2. (color online) Schematic diagram showing orienta¬ 
tional and vibrational degrees of freedom. The top right inset 
shows the potential energy landscape with (dashed red line) 
and without (solid black line) a cavity. The other insets show 
(top) the orientational degree of freedom and (bottom-right) 
the vibrational degree of freedom illustrating how the cou¬ 
pling to light decreases the displacement between the ground 
and excited manifolds. Shifts due to coupling to the cavity 
are exaggerated for clarity in all three insets. 


the matter-light coupling coupling arises due to the exis¬ 
tence of the counter-rotating terms in the original Hamil¬ 
tonian: If these terms were neglected then the energy 
in the ground state sector can be trivially found as the 
ground state would correspond to the empty state, and 
its energy would therefore not involve the matter-light 
coupling strength at all [s^. The presence of counter¬ 
rotating terms means that the ground state sector also 
involves an admixture of all even parity sectors, and the 
degree of admixture depends on the effective matter-light 
coupling. 

As discussed below, while exact solutions are possi¬ 
ble in some limiting cases of the orientational problem, 
these are not generally possible at finite temperature, 
nor for the vibrational problem. This is because ther¬ 
mal or quantum fluctuations of the auxiliary degrees of 
freedom break translational symmetry, preventing sim¬ 
ple exact diagonalization. As such, we proceed using the 
Schrieffer-Wolff formalism [4l|, which allows us to con¬ 
sider perturbatively the effects of these counter-rotating 
terms, and how they modify the energy landscape seen 
by the auxiliary orientational or vibrational degrees of 
freedom. For completeness, appendix |B] provides a brief 
summary of the Schrieffer-Wolff formalism. The essen¬ 
tial point is to separate H = Hg + Hi, where Hi are the 
terms treated perturbatively. At leading order this gives 
an effective Hamiltonian: 


H 




G,Hi 


G : [G,Ho] = iHi. (9) 


Taking the counter-rotating terms as Hi, the pertur¬ 
bation theory is controlled by the small parameter 
5 k,n/(wk + €n), which is indeed small for the physical 


parameters we consider [d^. We are thus in a regime 
where the counter-rotating terms cannot be ignored, but 
where they can be included perturbatively. In the fol¬ 
lowing sections we apply this approach in turn to the 
orientational and vibrational degrees of freedom, and see 
how the effective matter-light coupling can be derived 
self-consistently. 


A. Orientational degrees of freedom 

As discussed above, we consider first the classical 
orientational degrees of freedom subject to a pin¬ 
ning potential A„(0„). If all molecules are identical, 
An(0n) = A(0„), then one can find the zero tempera¬ 
ture ground state by choosing On = 0. For the ground 
state, all that is required is to find the quantum ground 
state of Eq. ([7]) as a function of 0 and then minimize over 
0. Since Eq. © is translationally invariant in the case 
On = 0, it is possible to find the exact ground state in 
the bosonic approximation by Fourier transforming. The 
bosonic approximation assumes the occupation of each 
excited molecule is small, a result that is valid unless 
5 k ^ Wk, e [H, 11^. However, at finite temperature, even 
for identical molecules, it is crucial to allow independent 
fluctuations of each assuming = 5 massively un¬ 
derestimates the entropy at finite temperature. At zero 
temperature, one may compare the exact solution to the 
Schrieffer-Wolff perturbative expansion used below, and 
one finds that these indeed match to leading order. 

For this rotational case, the form of G required 
in Eq. can be found trivially, and the resulting 

Hamiltonian can most conveniently be written as H = 

X^k ^k'0k'0k + with the molecular Hamiltonian 

in the form: 

% ^ )j(0,RWA) _ gk,nCOS^(6>„) 

4 ^ Wk + e„ ’ 

k 

where is the bare molecular Hamiltonian, in¬ 

cluding the RWA coupling to light and including the pin¬ 
ning term A„(5„). 

In order to consider the thermal distribution of 0„, 
we must specify the orientational potential A„(5„). We 
consider a form 

A„(5„) = Asin2(^^^-^^, (10) 

which tries to pin the molecules at angle 0°, relative to 
the cavity electric field, with strength A. We may thus 
consider both the free orientation case, A = 0, and the 
pinned case simultaneously. In what follows we will as¬ 
sume that the pinning angles have a uniform distribution 
such as would be found in a polymer matrix. This treat¬ 
ment, however, ignores effects which would be important 
in systems such as organic crystals in which the pinning 

























6 


angle has a fixed direction. The energy landscape for Defining = ((cos^(0„))), this ratio can be cal- 

each angle 9n, given the pinning angle 0° is then: culated as 


£;(0„|0°) = Asin2 


/ On - gl COS^ (On) 

V 2 J ^ Wk + e 


( 11 ) 


^ ^ J: _ ±_ 

g'^ 27r/3 dK^ 


d0°ln[Z(A,Ko,0°)], 


(15) 


where for simplicity we have assumed that each molecule 
has the same values of e, and the only disorder is in 
the pinning angles. In the following we will define the 
quantity: 


where the partition function is 


2{X,Ko,e°) 


= dO exp 


PKq cos^ 9 — p\ sin^ 




K„sY, 

k 


0k 


(wk + e) 


( 12 ) 


The quantity Kq characterizes the self-consistent energy 
favoring alignment of molecules. Replacing the summa¬ 
tion by an integral, and inserting the explicit forms of gk 
and Dk written above we have that 


As noted above l3Ko <C 1, and so we may Taylor expand 
in this small parameter to get a closed form for ^eff- As¬ 
suming that the pinning angle distribution = l/27r 

is uniform we obtain the simple result for the coupling 


1 r ^0 / J2(/3A/2)^ Y 

g2 2 [ ^ 4 V wmy. 


(16) 


/eC kdk 1 

Ic Jo ( 27 r)wfc(wfe-|-e)’ 


(13) 


where Ic is again the cavity length, C the combination 
defined following Eq. (ED, and A = 27r/aBohr is a cutoff 
reflecting the breakdown of the dipole approximation. To 
evaluate such integrals it is useful to write the dispersion 
in the form uj^ = (wg +c^k'^) + A(Njn/V which allows us 
to find the exact result 


with In{z} the imaginary Bessel function 

In[z) = jd9cos{n9)e^^°^^. (17) 

At this point we have made no assumption about the 
pinning strength A, and the value of the effective coupling 
is controlled by the combination /3A. If the pinning is 
strong /?A » 1 then we find the asymptotic form 


Ko = 


M 

2'kIcC^ 


In 


/ e -I- cDa A 
V e + Wo / 


(14) 




(18) 


A notable feature of Eq. da, as anticipated above, 
is that this quantity does not simply increase as one in¬ 
creases the polariton splitting by varying the density of 
emitters, Nm/V . There is a sub-leading dependence on 
the molecule density, via the renormalization Wk Wk 
given in Eq. O- However this effect is only significant 
in the deep strong coupling limit Physically 

this lack of scaling with Nm is because this “molecular 
adaptation energy” depends on the shift seen for each 
molecule. Inserting typical experimental values for an 
organic system Njn/V = 4.2 • 10^^m“^, e = wg = 2.1eV, 
Ic = 145nm, / = 0.5, OBohr = 2nm, C = 4.3eV^nm^ 
which correspond to the values extracted from Ref. [^ . 
one finds that 

Ko = 6.3 • 10"VeV, 

which is much smaller than ksT at room temperature. 

As noted earlier, the effective matter-light coupling 
strength is given by g^^k = 0k((‘^o®^(^")))- 
temperature, this corresponds to minimizing the energy, 
leading to A„ sin(0„ — 0°) = —2Arg sin(20„). However, 
since the value of Kq given above is such that Kq <C ksT, 
it is crucial to consider finite temperatures. The small¬ 
ness of the ratio Ko/kgT will also allow us to make fur¬ 
ther perturbative expansions in the following. 


This no longer depends on temperature as the strong 
pinning limit means entropy become unimportant. Thus, 
in the A —> oo limit the coupling takes on the isotropic 
value of 1/2, corresponding to the uniform distribution 
of angles 0„ = 0°. The effective coupling increases as the 
pinning A decreases. In the limit of vanishing pinning 
A —>■ 0, the imaginary Bessel function /2(0) vanishes and 
so we have: 

f = + (19) 

Both Eq. (fT^ and Eq. (fl^ indicate that as long as 
fdKo “C 1, the modification and temperature dependence 
of g^g is very small. We note that while in the organic 
systems which we focus on here Kq is relatively small, in 
systems which have very small mode volumes this param¬ 
eter could be engineered to be much larger and hence the 
coupling strength renormalization would be much more 
pronounced, i.e., since there is no scaling with the 
crucial feature to see a strong renormalization is to min¬ 
imize the mode volume in absolute terms, and not the 
mode volume per molecule. This sug gest s evanescently 
confined radiation modes in plasmonic |45| or phonon po¬ 
lariton ( 4 ^ systems may be a promising venue to explore 
this physics. 
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B. Vibrational degrees of freedom 


In the vibrational case, even without disorder, an exact 
solution of Eq. ([5]) via Fourier transformation is no longer 
possible, because the vibrational degrees of freedom are 
modeled as quantum degrees of freedom with their own 
quantum dynamics. As such, they break translational in¬ 
variance — i.e. localized vibrational excitations can scat¬ 
ter between different polariton momentum states. Thus, 
once again we must use the Schrieffer-Wolff formalism. If 
we start from Eq. ([5]), solving the equation [G, Hq] = iHi 
is now more challenging than for the rotational case, as 
Hq involves terms that couple the electronic state to the 
vibrational quantum state. The equation can however be 
solved in the form of a power series, 




-WKr. 


k,n,j 


(Wk + 


- H.c.) , (20) 


where the operators Oj are defined by the recursion re¬ 
lation 


m 






dj-iy^^{b 

n 



with the base case Oq = 1. This expansion then allows 
one to write out the effective Hamiltonian in the same 
form as above, but with the molecular Hamiltonian, 


hn 


_ 1^(0.RWA) 





^3 


( 21 ) 


with the bare molecular Hamiltonian, includ¬ 

ing the RWA coupling to light, and vibrational terms of 
Eq. (0). 

This expression can be considered as a multinomial 
power series in the quantities Hm/(wk + e„) <C 1 for 
each vibrational mode m. For typical parameters, such 
quantities are small, and so we may truncate at first or¬ 
der, i.e. keep terms up to j = 1. Beyond j = 1, the 
expression becomes considerably more complicated, as 
cross terms between different vibrational modes appear. 
Up to j = 1, we find an effective molecular Hamiltonian 
which we write out in full (neglecting constant terms): 


f Cn + Kq 


E 




cr„ -b H.C. 


(blb^ + ^ibl + Ly 

m ' 

'JSm{bm + ) ■ ( 22 ) 




I - a. 


The j = 0 term gave an energy shift of two-level systems 
with the same form Kq found previously. The j = 1 term 
is on the last line, and describes a shift to the vibrational 
modes. The coefficients for these terms are defined by a 
generalization of that used in Eq. (fT^ for the rotational 
case above. 


K,y) = Y. 


5k 


(wk + e)^+i ’ 


where we may find an analytic form for the resulting 
integral 

I 1 

(e -|- wo)'l (e + wa)-* 


^ 2'KlcC^j 


We may once again note that none of the terms Kj (e) are 
proportional to the number of molecules Nm'- vacuum- 
state molecular adaptation is not collectively enhanced. 

From the form of Eq. (1221) it is clear that the term in Ki 
describes a reduction in the offset between the vibrational 
ground state in the two electronic states, as the virtual 
excitations admix the excited electronic state configura¬ 
tion into the ground state. This can be viewed as a reduc¬ 
tion of the Huang-Rhys parameter, S'm.eff = S'm)! —2iGi]. 
This point is illustrated in Fig. [2l As noted earlier, the 
dependence of on the vibrational degrees of freedom 
is more complicated: one must calculate the overlap be¬ 
tween the vibrational states of the ground and excited 
electronic state manifolds. The reason that the vibra¬ 
tional states differ in these manifolds is the existence of 
the terms {by+b^)(j^^ which correspond to displacement 
of the vibrational coordinate dependent on the electronic 
state of the molecule. As such, a reduction of Sm,eS would 
lead to an enhanced matter-light coupling. However, us¬ 
ing the same typical experimental values as before we find 
that this dimensionless shift has a value of approximately 
Ki = 2.9 • I0“®. Thus, once again one may conclude 
that the characteristic scale of any vibrational molecular 
adaptation (determined by Ki) is negligibly small. 


C. Other aspects of molecular state 

The discussion so far has focused on two specific 
microscopic mechanisms which might have led to self- 
consistent adaptation of the molecules so as to enhance 
their coupling to light. Here we note those aspects of the 
above results which can be easily generalized to other mi¬ 
croscopic mechanisms. Examples of such other potential 
mechanisms include solvation of the molecule, molecu¬ 
lar configuration, and charge transfer state. In some 
cases, these will require detailed modeling of the spe¬ 
cific process, particularly for degrees of freedom with en¬ 
ergy scales larger than temperature, where quantum ef¬ 
fects become important, as in the example of vibrational 
modes above. However, the basic idea can be illustrated 
in the simplest case, where some aspect of configuration 
can be parametrized by a classical variable. 



















In the case where classical parametrization applies, the 
description is very similar to the discussion of orienta¬ 
tional degrees of freedom: One considers a classical vari¬ 
able Xn which parametrizes some aspect the state for 
molecule n. Associated with this will be an energy func¬ 
tion Kn(xn), and, for self-consistent adaptation to oc¬ 
cur, the matter-light coupling must depend on this vari¬ 
able as g ^ g4>n(xn)- The same perturbative analysis 
as discussed above will then lead to the effective energy 
function: En{xn) = A„(a:„) — A"o(^n(a:n)^- This in turn 
means that the self-consistent effective matter-light cou¬ 
pling will take the form: 


_ 1 ^ In 

gl PdKo 


dxr,e 




Since this involves the same sum Kq as defined in 
Eq. (fT^ . the same absence of scaling with number of 
molecules occurs i.e. it is the single-molecular, rather 
than collective, coupling which is important. Moreover, 
since PKq <C 1, then: 


9l,eS - 9l{{(l>n{Xnf))o + 0{/3Ko) 

where ((.. .))o indicates thermal averaging with the bare 
energy function An{xn), i.e. any such vacuum-state self- 
consistent adaptation of molecules is suppressed by the 
small parameter /3A'o- 


III. VIBRATIONAL DRESSING OF THE 
EXCITON-POLARITON SPECTRUM 


In the previous section we have seen that while self- 
consistent molecular adaptation due to matter-light cou¬ 
pling is possible, neither the strength nor the dependence 
upon molecular concentration are consistent with the re¬ 
sults reported in Ref. [2^. In this section we show how 
directly calculating the absorption spectrum from the 
model of vibrationally dressed polaritons introduced in 
the previous section can lead to a rather different mecha¬ 
nism which could however be responsible for the features 
observed. This mechanism arises from the combination 
of vibrational dressing of the spectrum and the effects 
of disorder. We focus in particular on the peak in the 
spectrum near the bare excitonic resonance. We will see 
that the shape of this feature depends strongly on the 
vibrational state of the molecules. 

This section is divided into two subsections, in sec¬ 
tion IIII Al we first discuss how vibrational excitations 
should be included in the calculation of the disordered 
polariton spectrum. This follows the method outlined in 
appendix |B1 and so all that is required is to calculate the 
excitonic self energy. We then discuss the resulting form 
of the spectrum in section UlI Bl 


A. Self energy of vibrationally dressed excitons 


As discussed in appendix lAl the absorption, emission 
and transmission spectra can all be found in terms of 
the photon retarded Green’s function. In this approach, 
all the properties of the molecules (i.e. inhomogeneous 
broadening, vibrational dressing etc.) are incorporated 
via the excitonic self energy I^k,xxiE). We must there¬ 
fore calculate this quantity, defined in Eq. (IA1IA2II . for 
the vibrationally dressed Hamiltonian, Eq. ([8]). To do 
this, it is useful to label the eigenstates as Ipt);!'?)-): 
corresponding to the vibrational eigenstates in the ex¬ 
cited electronic state manifold and the ground electronic 
state manifold. We denote the energies of these states 
as Af, and and we introduce the overlap matrix el¬ 
ement: apg = This matrix element de¬ 

scribes the extent to which the state with q vibrational 
excitations in the electronic ground manifold overlaps 
with the state in the excited electronic manifold with p 
excitations. In order to find these eigenvalues and eigen¬ 
functions we must diagonalize the vibrational problem, 

Hvih = -I- cr^ in the electronic ground 

and excited states. As discussed in Sec. IIIBl the renor¬ 
malization of these parameters due to virtual pair cre¬ 
ation is very small, and so we neglect it in the following. 

In terms of these quantities we may write the self en¬ 
ergy, including inhomogeneous broadening of excitonic 
energies in the form: 


vRWA 


(^) 


9ti^m 

Z 


d£h{e) 


0!.pq 



ly + ij + {Eg - e],) 
(23) 


where h{e) is again the distribution of exciton energies, 
as in Sec. IE and we have again ignored the subleading 
effects of the exciton energy distribution on the coupling 
strength, by using g^Nj^ = ffk n- should be noted 
that the energies Ej,,E^ depend (linearly) on the energy 
e, but that the matrix elements apg are not dependent 
on this energy scale. 

In the following, we model the inhomogeneous broad¬ 
ening of excitons by using the distribution h(e) oc 
0(e)e*^'^“'’°^ , i.e. a truncated Gaussian distribution of 

excitonic energies. The truncation has little effect, but is 
formally required, as negative energy states cannot phys¬ 
ically exist, and cause problems in the ultrastrong cou¬ 
pling limit [13,113 ■ 


B. Dependence of absorption spectrum on 
vibrational state 

The exciton spectrum in the absence of vibrational 
dressing was shown previously in Fig. [1] Before exploring 
the effect of vibrational modes, it is helpful to summarize 
how the residual exciton peak arises. Mathematically, 
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the excitonic feature in the polaritonic spectrum can be 
understood directly from the form of Eq. (0 and the 
definition of the absorption spectrum. The peak at the 
exciton frequency occurs because the imaginary part of 
the retarded Green’s function has a numerator involving 
the imaginary part of the self energy, and 

this self energy has a peak at the excitonic energy. How¬ 
ever, the weight of the residual excitonic peak reduces as 
the matter-light coupling increases, because the excitonic 
self energy also appears, squared, in the denominator. 
It is also important to note that this residual excitonic 
peak does not appear in the transmission spectrum; since 
T\^{v) oc \G^, there is no term in the numerator 
of T\^{v) arising from the exciton self energy. Thus, the 
residual exictonic peak is a feature of the absorption (and 
reflection), but not the transmission spectrum. 



v[eV] 



FIG. 3. (color online) The evolution of absorption spectra 
with the strength of gy/N for a system with a single vibra¬ 
tional mode with S = 0.06, Q. = O.OSeV. The insets show the 
equivalent absorption spectra for bare molecules. The other 
parameters are as in Fig. [T] Panel (a) shows the spectrum 
for a perfect cavity while panel (b) includes the effects of 
cavity losses and non-radiative excitonic decay exactly as in 
Fig. mb). 

Figure [3] shows an equivalent set of spectra to Fig. [U 
but with vibrational dressing. For comparison in the in¬ 
sets we show the absorption spectra of a bare molecule 
Q'[Sk,xx(i^)], i-e. the spectra without a cavity. Panel (a) 
shows the simple case where cavity losses are ignored, 
while panel (b) shows the more experimentally relevant 
case where these effects are included (discussed in ap¬ 
pendix In the following we use the notation p — q 
which denotes transitions in the molecule from the state 
with p vibrational excitations in the electronic ground 
state to the state with q vibrational excitations in the 
electronic excited state. The presence of the vibrational 
modes causes dramatic changes to the residual excitonic 
peak in the polariton spectrum not seen in the bare ex¬ 
citonic absorption. For the bare molecule the spectral 


weight associated with transition with the “zero phonon 
line”, i.e. the transition denoted 0 — 0 in the notation 
introduced above is completely dominant, only a small 
amount of weight visible in the sideband which corre¬ 
sponds to 0 — 1 transition. When the molecule is placed 
inside a cavity, the vibrational sidebands corresponding 
to the 1 — 0 and 0 — 1 transitions become much more 
prominent. Physically this is because the spectral weight 
that had been associated with the 0-0 transition in the 
bare molecular spectrum has been moved into the polari¬ 
ton peaks of the spectrum. The mathematical form of the 
Green’s function makes clear that formation of the polari¬ 
ton spectral feature is predominantly at the expense of 
whatever feature dominates the excitonic emission spec¬ 
trum; here this is the 0-0 feature. The excitonic feature 
corresponds to the “left-over” spectral weight associated 
with the subradiant states. Thus, the vibrational side¬ 
bands have been “excavated” by removing the dominant 
feature from the zero-phonon line. Including the the ef¬ 
fects of cavity losses and non-radiative excitonic decays, 
as can be seen in Fig. [3] (b), washes out this complex 
sideband structure of the central peaks and the spectrum 
as a function of coupling strength looks very similar to 
that obtained without coupling to vibrational modes, as 
in Fig. [T] However as we discuss below, there are still 
effects which can be observed which are a direct conse¬ 
quence of the vibrational structure. 

Fig. m illustrates the evolution of spectra with tem¬ 
perature. Panel (a) shows that without coupling to vi¬ 
brational modes there is no notable temperature depen¬ 
dence. In the presence of the vibrational dressing, a 
strong temperature dependence appears. At higher tem¬ 
peratures, there is a greater thermal occupation of the 
vibrational modes hence the spectral weight under the 
vibrational peaks rises. While this has a small effect in 
the bare molecular spectrum (where the 0-0 transition 
dwarfs all other features), see inset in (b), it is very pro¬ 
nounced in the polariton spectrum and is even visible in 
the presence of large cavity losses as in Fig. HKc). 

The figures so far have shown results where disorder is 
relatively small, and so vibronic replicas can be clearly 
observed for the good cavity, but merge for the bad cav¬ 
ity limit. To show that small disorder is not required 
for the strong temperature dependence to occur Fig. [5] 
shows the effects of large disorder (i.e. inhomogeneous 
broadening). Just as seen for the homogeneous broaden¬ 
ing in Fig. [UJc), a temperature dependence of the resid¬ 
ual excitonic peak is still visible. In this figure, we have 
also included a more complicated vibrational spectrum, 
involving two vibrational modes. Such a system will ex¬ 
hibit behavior similar to that of a single mode but with 
a large vibrational coupling . 


IV. DISCUSSION 

In this paper we have presented two microscopic mod¬ 
els which could in principle describe self-consistent molec- 
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FIG. 4. (color online) Evolution of absorption spectra with 
the temperature (values shown on the colorscale) for: (a) A 
system with no coupling to vibrational modes, S = 0, (b) Cou¬ 
pling to a single vibrational mode with S = 0.02, Q — O.OSeV. 
The inset shows the bare molecular spectrum, (c) Including 
the effects of cavity losses and non-radiative excitonic decay. 
All panels are for gy/N — 0.3eV, and other parameters as in 
previous figures. 



FIG. 5. (color online) Temperature dependence of absorption 
spectrum for a system with strong disorder a — 0.025eI4. The 
main panel shows the optical spectrum of the strongly coupled 
system, the inset shows the spectrum of the bare excitons. In 
this case, we include two vibrational modes: = 0.2, Hi = 

O.OSeV, and S 2 ~ 0.3, H 2 = 0.04eV. Other parameters as in 
Fig.H 


ular adaptation so as to maximize the vacuum-state cou¬ 
pling to light. In both cases, the crucial feature of the 
model is the counter-rotating terms in the matter-light 
coupling. These allow virtual fluctuations in the ground 
state, that lower the ground state energy depending on 
the configuration of the molecules. This energy gain is 
the only energy gain that can be relevant in the linear 
response regime — i.e. the question of whether the ex¬ 
cited states would have lower energy is not of relevance 
while the system is only weakly pumped. We found that 
while such a mechanism for molecular adaptation does 
exist, it does not show any collective enhancement, in 
contrast to the polariton splitting, and does not there¬ 
fore lead to significant molecular adaptation, even when 
the polariton splitting gky/Nm ~ w, e. 

The appearance of a residual exciton peak in the po¬ 
lariton spectrum would be affected by any such self- 
consistent molecular adaptation if its scale were suffi¬ 
cient. However, for relevant parameters, such effects are 
dwarfed by a far more dramatic effect, of vibrational 
dressing of the residual exciton peak. This leads to a 
pronounced temperature dependence of the feature near 
the exciton energy in the absorption spectrum. 

While the molecular adaptation energy scale is not col¬ 
lectively enhanced, the basic underlying physics could po¬ 
tentially be relevant in single molecule strong coupling, 
e.g. with plasmonic resonances [i^. In such cases, rather 
than having many molecules in a large mode volume, the 
mode volume is reduced so that strong, or even ultra¬ 
strong coupling occurs at the single molecule level, so 
that both the polariton energy and the molecular adap¬ 
tation energy become large. In such a case, one may 
hope to see either reorientation, or renormalization of 
the Huang-Rhys parameter due to strong coupling. 

Another intriguing direction for future research is to 
consider how the physics discussed in this manuscript 
interacts with the physics of polariton condensation and 
lasing [5^ . Polariton condensation has been seen in both 
inorganic [5l|, [111 and organic fTl - [l8l| systems. In addi¬ 
tion, condensation of photons has been seen for weakly 
coupled sterns of organic molecules [5l|. Theoretical 
work [3; has begun to address some of the pecu¬ 

liarities of the organic polariton system, including effects 
of disorder and of vibrational modes. However, features 
as seen in this paper, resulting from the interplay of these 
may lead to further exotic behavior in the high density 
condensed phase. 
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Appendix A: Absorption, Transmission and 
Reflection spectrnm of exciton-polariton system 

In this appendix we summarize the calculation of the 
absorption, transmission and reflection spectra. Some 


subtleties arise because we wish to calculate the spec¬ 
trum of a model with ultrastrong coupling, i.e. without 
making the rotating wave approximation. Such results 
were first calculated by Ciuti and Carusotto [i^, here 
we present a synopsis of these results, as well as a “dic¬ 
tionary” to translate the results of that paper into the 
language of Green’s functions. We begin by defining 
the retarded Green’s function [^. Because we consider 
both CO- and counter-rotating terms, we must consider 
both normal and anomalous Green’s functions, i.e. we 
must include number non-conserving terms which appear 
for ultrastrong coupling, and thus we consider a matrix 
Green’s function; 

r^(f f'] = -i ( ([V’lkW>V’k(^')]) \ 


I 

In terms of the Bogoliubov transformed operators, i.e. the operators appearing in Eq. m the inverse Green’s 
function takes the form 




fi^ + ik{v) - Wk + +ih(v) -I- Ek,xx(i^) 


where xxi^) energy for a photon of in-plane momentum k, arising from the excitonic response (discussed 

further below), and k{v) is the loss rate. 


Here, following we have used a frequency depen¬ 
dent complex loss rate k(y). Frequency dependence is 
required for physical consistency in the case of ultra¬ 
strong coupling — Markovian loss and ultrastrong cou¬ 
pling would predict a perpetual light source. Frequency 
dependent loss requires, via the Kramers-Kronig relation, 
a corresponding Lamb shift, which is incorporated into 
the imaginary part of k{y). Both S and k are written 
for the Bogoliubov transformed operators, and so both 
these terms incorporate a pre-factor Wk/wk to account 
for the Bogoliubov transformation of the combination 

(4+^lk)- 

The specific self energy required, Ej, {v ), corresponds 
to correlation functions of the excitonic operators. In 
the absence of strong (i.e. beyond RWA) excitonic damp¬ 
ing (see for the more general case), this self energy 
can however be related to results in the rotating wave 
approximation by: 


self energy, 


vRWA 


{iUJm) 



^^{p\^n{T)cr^(Q)\p)e (A2) 


and replacing the Matsubara frequency by iuJm v + 
i0+. The sum over p appearing here is over all states 
of the excitonic system, and Z is the partition function. 
For the “vacuum” state we consider — i.e. in the absence 
of strong pumping — these are the bare exciton states, 
including the quantum states of any auxiliary degrees of 
freedom. 

Using the input-output formalism [fl^ adapted to the 
ultrastrong coupling regime [1^, one can write a fre¬ 
quency dependent scattering matrix relating input and 
output fields at the left and right sides of the cavity, 


= A (eK*M + RT_(-)]*). (Ai) 

where the expression E^'^^ is the “standard” self energy 
that would appear in the rotating wave approximation, 
depending on the correlation of a~ operators. These 
can most easily be found by analytic continuation from 
imaginary time to real time, starting from the Matsubara 


1 - »d>k,L(t^) -*\/0k,L(jA)(/)k,fl:(!^)\ 

-i\/ I - i(l)k,R{i') ) 

where xxi’^) fl(^) the real 

parts of the loss rates arising from the left and right mir¬ 
rors and the quantity relates to the matrix re¬ 
tarded Green’s function as = l'^G^{i')I where 

= (l !)■ This structure means the Bogoliubov trans¬ 
formation corresponds to = (oJk/^k)G^^^(i') 
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and the Bogoliubov transformed Green’s function takes 
the form: 

^^ -■ (A3) 

’ 1/2 _ ( 2,2 _|_ 

One can then find the transmission Tk(i^), reflection 
R\^(y) and absorption Ay^{h') coefficients by consider¬ 
ing the modulus square of various coefficients. Clearly 

is independent of which 
direction light is incident from, while the absorption co¬ 
efficient ^k,cr=L,ij(' 2 ) takes the form 


AkA^) = -k'^{v) [23 [G^^^^iv)] + k{v)\G^^^^{v) 


. 

with k{v) = {kl{v) + k,r{v))/ 2 . The prefactor in this 
expression shows the obvious dependence on the trans¬ 
missivity of the input mirrors. 

In order to separate mirror transmissivity dependent 
features from the “intrinsic” properties of the ultra¬ 
strong coupling we will consider below the two quantities 
^k(^) = |G« as being proportional to the trans¬ 
mission, and ak(i^) = —23[Gk as controlling the 

absorption in the limit of a good cavity, i.e. k{v) —^ 0 . 
The quantity a]s_{v) differs from the full absorption as it 
neglects interference effects at the input mirror. 



V [eV] 


FIG. 6 . (color online) Absorption spectrum, as defined by the 
full expression, Eq. (TAil) . plotted for = 0.25, S 2 = 0.35, 
fli = O.OOeV, 0,2 = 0.05eV, ksT — 0.035eV and various cav¬ 
ity loss rates, n. Other parameters are as in Fig. As dis¬ 
cussed in the text, in order that the polaritonic peaks have a 
finite width, one must include the effects of excitonic absorp¬ 
tion or non-mirror cavity losses. For this figure, we include 
an excitonic linewidth 7 = 10“‘*eV. 


For comparison to the results in Ref. [ 2 ^ where silver 
mirrors were used, we present in Figures [H [31 and 0] the 
effects of large cavity linewidth. Figure [ 6 ] also shows how 
the full absorption spectrum given by Eq. (IA4I) evolves 
with varying linewidth. In calculating these spectra for 
large linewidth, an issue arises regarding the form of the 
absorption spectrum: The equations written above as¬ 
sume that the only photon loss is due to escape through 
the mirrors. This means that the k{u) appearing explic¬ 
itly in Eq. (IA4D (describing interference effects from the 


mirror) is the same as the K(y) appearing in the denomi¬ 
nator of the photon Green’s function, Eq. (ESI), and one 
may check that for any frequency where S[Ek,a;a:(j^)] is 
small, this causes a near cancellation between the two 
contributions. For the Gaussian exciton density of states 
h{e) used in this paper, this cancellation almost com¬ 
pletely suppresses the polariton peaks. Such a cancel¬ 
lation in the absorption spectrum can be expected on 
physical grounds: if the only loss channel for photons 
is the mirrors, then there is no absorption. All photons 
that enter eventually leave. In a real device there are 
other photon loss sources (absorbers, scattering by sur¬ 
face roughness). Similarly, in a real device, the excitons 
have a non-zero rate of non-radiative decay. The effect 
of this is included by retaining a non-zero value of 7 in 
the denominator of the self-energy, Eq. (ESI). This leads 
to Lorentzian tails of 5[Ek,a:a;(i^)], giving a finite weight 
to the polariton peak in the absorption spectrum; such 
an effect was included by Houdre et al. [S^. We follow 
this approach in plotting Eig. [HI 


Appendix B: Schrieffer-WolfF transformation 


In section El we make use of the Schrieffer-WolfF ap¬ 
proximation; for completeness we provide here a brief 
explanation of this formalism. The approach is based on 
dividing the Hamiltonian into two parts, H = Hq + Hi, 
where the term Hi takes one between different “sectors”. 
In our case, these sectors correspond to different numbers 
of polaritons — i.e. Hi is the “counter-rotating” part 
of the Hamiltonian which simultaneous creates a photon 
and excites a molecule. 

The aim of the Schrieffer-Wolff formulation is to make 
a unitary transformation H = such that the 

transformed Hamiltonian no longer has any coupling be¬ 
tween sectors. Physically, this corresponds to eliminat¬ 
ing the effect of virtual pair creation and destruction, 
and deriving how such virtual processes renormalize the 
Hamiltonian within a given sector. 

If the Hamiltonian Hi can be treated perturbatively 
by replacing Hi ^ rjHi with 77 a small parameter, then 
one can consider a series solution G = . In 

order to make the first-order terms in 77 vanish, one must 
choose [G(i),Ho] = iHi. This then leads (setting 77 = 1) 
to the expression: 


H = Ho + - 


2 1 


Gil), Hi 


H.o.t 


(HI) 


where the higher order terms involve all Stop¬ 

ping at leading order gives the expression in Eq. m , cor¬ 
responding to the leading order effects of virtual pair cre¬ 
ation and annihilation. To solve [G^Hq] = iHi in prac¬ 
tice is straightforward if one knows the eigenspectrum of 
Hq = En'^\'n-){n\, which then allows one to write 


G=iY: 

n,m 


\n){n\Hi\m) {m\ 
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